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Abstract 

IFS Fractals - the attractors of Iterated Function Systems in Euclidean space - have 
motivated plenty of investigation to date, due to their simplicity and applicability in 
various fields, such as the modeling of plants in computer graphics, and the design of 
fractal antennas. The statement and resolution of the Fractal-Line Intersection Prob- 
lem is imperative for a more efficient treatment of applications. This paper intends 
to take further steps towards this resolution, building on the literature. We provide a 
verifiable condition guaranteeing intersection with any line passing through the con- 
vex hull of an IFS fractal, shown in M. d for hyperplanes. The condition also implies a 
constructive algorithm for finding the points of intersection. We give several numerical 
conditions that guarantee an infinite number of approximate intersections if there is at 
least one. In our effort to quantify the intersection, we provide an explicit formula for 
the well-known invariant measure of IFS.* 

MSC class: 28A80 (primary); 37F99, 52A35 (secondary) 
Keywords: fractals, attractors, IFS, invariant measure, transversal 

Falconer [7] provides an extensive overview of Fractal Geometry, summarizing the work of 
several researchers. In particular, he examines the behaviour of Euclidean attractors under 
projections to hyperplanes. When the projection is carried out in M 2 or IR 3 onto a line or 
plane, we can think of the resulting set as the "shadow" of the fractal, which may be analyzed 
for its own dimension. Furthermore, we may ask how many fractal points are projected to a 
certain point on the line, or if there are any at all, in essence trying to uncover the distribution 
of such a projection. But what is the fundamental underlying question here? Upon further 
thought, we may realize that the projection problem breaks down into two main questions, 
since the directional ray of projection can be thought of as an intersecting line. 

*A preprint of the article has been submitted for consideration in Fractals © 2012 World Scientific 
Publishing Company http://www.worldscientific.com/worldscinet/fractals 
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Therefore we must examine the following questions: 

(1) The Fractal-Line Intersection Problem: Given a line and an IFS fractal in the plane, do 
they intersect? 

(2) If they intersect, how many points of intersection are there? How about within some 
e > neighbourhood of the line? 

Our paper intends to provide resolutions to these questions, while demonstrating their rel- 
evance from an applicational point of view. We show that for certain broad classes of IFS 
fractals - hyperdense, or specifically chain fractals - the shadow is always filled in, no matter 
where the light shines from, proven in general in IR d . Furthermore, we show that the segment 
shadow in IR 2 receives an infinite number of projected points, in any e > subinterval of the 
segment. These properties make some of these potentially disconnected fractals ideal for 2D 
fractal antenna design or as light-absorbing tree crowns (suggesting that tree crowns could 
be considered 3D fractal antennas). 

Recently Mendivil and Taylor [20] have approached these questions from a projectional point 
of view. Defining a certain parametrized class of planar IFS fractals, they wish to guarantee 
that the shadow in all directions is a segment. They prove that this holds for some domain 
of parameters. In other words, for this specific class of planar fractals, any line or ray of light 
that intersects the convex hull, also intersects the fractal, thereby contributing to its shadow. 
We examine the problem further for the broadest possible class of attractors in IR d called 
hyperdense fractals, and we also introduce the verifiable subclass of chain fractals. 

The Fractal-Line Intersection Problem is relevant to a number of applications, among which 
we mention ray tracing in computer science, the design of fractal antennas in engineering, 
and the study of tree crown density for light absorption in botany and forestry. In computer 
science, ray tracing involves the shading of an object in virtual space, which is detailed in 
Hart and DeFanti [10] in regards to 3D IFS fractals, as well as numerous other publications. 
Fractal antennas are flat metal antennas with an IFS fractal layout that must be optimized 
for the amount of material used vs. the efficiency of signal reception. These antennas were 
introduced by Cohen et al. [2, 11]. The study of tree crowns for light absorption is a vast 
field, and we recommend the work of Lindenmeyer et al. [15] and Zeide [21]. Last but not 
least, we emphasize the pioneering inspirational work of Mandelbrot [16]. 

Further research into the projection of fractals has been carried out by Besicovitch [1] and 
Federer [9] examining s-sets; Marstrand [17], Kaufman [14], and Mattila [18, 19] showing 
projection theorems for arbirary sets in IR d ; and Davies [3], Falconer et al. [6, 8], and Howroyd 
[12] giving results for box and packing dimensions. These efforts are all summarized in the 
expository book by Falconer [7], which also provides an introduction to IFS fractals. Most 
recently, Elekes et al. [4] investigated the reconstruction of geometric objects based on their 
intersection with test sets, emphasizing its relevance to tomography. 
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1 Preliminary Concepts 



The attractors of Iterated Function Systems - IFS Fractals - were first defined and exam- 
ined by Hutchinson [13], and may be the most elementary fractals possible. They are the 
attractors of a finite set of affine linear contraction mappings on IR d - the "function system" 
- which when combined and iterated to infinity, converges to an attracting limit set, the IFS 
fractal itself. We begin by defining the generating maps of an IFS Fractal, and go on to 
stating its existence and uniqueness. 

1.1 IFS Fractals 

Definition 1.1 Let an affine contraction mapping (briefly: contraction) T : R d -> Mr be 
defined for all z G IR d as T(z) := p + M(z — p) where p G IR d is the fixed point ofT, and the 
invertible M G R dxd is the factor ofT, with \\M\\ 2 < 1. 

When T is a similitude over the complex plane (d = 2), it may be written in the form 
T(z) := p + ip(z — p) where z,p G C and ip = \e dl G C, with A G (0, 1) and 9 G (— 7r, tt]. In 
higher dimensions, similitudes are M = XR with a unitary R. 

Definition 1.2 Let an affine contractive n-map iterated function system (briefly: IFS) be 
defined as a finite set of contractions, and denoted as T := {Ti, . . . ,T n }, n G N. We will 
further denote N := {1, ... ,n}, V :— {pi, . . . ,p n }, $ := {Mi, . . . , M n }. 

Definition 1.3 Let T = {Ti, . . . ,T n },n G N be an IFS. Define the Hutchinson operator H 
belonging to T as 

n 

H(S) := |J T k (S), T k (S) := {T k (z) : z G S}, for any S C C 
k=i 

and call H(S) the Hutchinson of the set S. 

Theorem 1.1 For any IFS with Hutchinson operator H, there exists a unique compact set 
Fctf such that H(F) = F. Furthermore, for any compact So C M d , the recursive iteration 
S n+ i := H(S n ) converges to F in the Hausdorff metric. 

Proof The proof follows from the Banach Fixed Point Theorem, once we show that H is 
contractive in the Hausdorff metric over compact sets [13]. □ 

Definition 1.4 Let the set F in the above theorem be called a fractal generated by an IFS 
with Hutchinson operator H (briefly: IFS fractal). Denote (Ti, . . . ,T n ) = (T) := F. 
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1.2 The Address Set 



The address set results from the iteration of the Hutchinson operator, and it is a way to 
label the location of each fractal point. Since we can start the iteration towards F with 
any compact set, we will choose the primary fixed point, which is any point in V of our 
preference. 

Definition 1.5 Let M'i :— Af x . . . x J\f be the index set to the j-th Cartesian power, and 
call this j the iteration level. Then define the address set as A := {0} U U^i-^ 7 - For any 
a E A denote its k-th coordinate as a(k), fceN. Let its dimension or length be denoted as 
\a\ E N so that a E 7V' a ' and let |0| := 0. Define the map with address a E A acting on any 
z ER d as the function composition T a (z) := T a ^ o . . . o T a ^ a ^(z). Let the identity map be 
T := Id. 

The above definition of the address set revises and unites the two standard definitions in use. 
Common ways include identifying addresses with decimal number representations with a cer- 
tain basis, and strings of letters usually for a small number of maps. Our formal language 
representation replaces the cumbersome decimal notation while accounting for finite ad- 
dresses, by using vectors of numbers preferable to the string formalism when \J\f \ > 24. 

Theorem 1.2 For any primary fixed point p EV we have 



We call this the address generation of F. 

Proof The proof follows from Theorem 1.1 with either of the initial sets {p} or V. □ 

Definition 1.6 Let the address adr(/) of a fractal point f E F with respect to a primary 
fixed point p EV be the shortest address a E A for which T a (p) = f . 

1.3 Hyperdense Fractals 

Definition 1.7 Let an IF S fractal in M. d be hyperdense if any hyperplane that intersects its 
convex hull, also intersects the Hutchinson of its convex hull. We will call an IFS fractal a 
chain fractal, if the Hutchinson of its convex hull is connected. 

Theorem 1.3 Chain fractals are hyperdense. 

Proof Denote the fractal F = (Ti, . . . ,T n ) and its convex hull Cf '■= Conv(F). First we 
see that for any TeTwe have T(Conv(S')) = Conv(T(S')) since T is affine. We also note 
that for any Si, . . . , Sn C M. d 



F = lim H N ({p}) = C\{T a (p) : a E A} = C\{T a (p k ) : a E A, p k E V} 



Conv I I Sf- I = Conv 
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since the convex combination of convex combinations, is a convex combination. Considering 
the fact that H(F) = F as well as the above ideas, we have that 

Conv(F) = Conv (\\J Conv(T fe (F)) j = Conv \ \J T fc (Conv(F)) J 

meaning that Cp = Conv (H(Cf)). Since F is compact so is Cf and thus H(Cf), since H 
is continuous. We now turn to showing that F is hyperdense. Let us take any hyperplane 
L C R d that intersects C F = Conv(H(C F )) in some point q = fih\ + (1 — /i)/i2, ^1,2 G 
H(Cp), H G [0,1]. If /x e {0,1} then g G H(Cp) so we are done. Otherwise L separates 
the space into two half spaces, with h% in one and h 2 in the other. Since F is a chain 
fractal, we know that H{Cp) is connected, thus it is path-connected, so there is a path 
7 C H{Cf) connecting h\ and /i2- Since hip are on separate sides of L, we must have that 
L PI 7 7^ 0. This can be shown by parametrizing 7 : [0, 1] -> IR d and writing the plane as 
L = {z G M. d : (a, 2) = 6} for some a G M d , 6 G R. Denoting f(t) = (a,j(t)} — b we have 
that /(0)/(l) < 0, so by Bolzano's theorem / must have a root £ £ (0, 1), implying that 

7(t ) e Ln 7 c LnH(C F ). 

Therefore by L fl 7 7^ and 7 C #(Cp) we have that L (1 H(C F ) ^ 0. □ 

Lemma 1.1 For a hyperdense fractal F = (Ti, . . . , T n ) and an?/ address a <E A, if a hyper- 
plane L intersects T a (Cp) then it also intersects T a {H{Cp)). 

Proof It is clear that since the M k G R dxd factors in T k are invertible, the inverses of the 
maps are T A T 1 (z) = p k + M^ 1 (z — p k ) . Thus T" 1 also exists, and it is also an affine mapping, 
so it takes the hyperplane L into another hyperplane V . Thus LC\T a (Cp) 7^ is equivalent 
to T~ X (L) H C F ^ 0, which by the hyperdensity of F implies that r a _1 (L) n H(C F ) ^ 0, 
and so L n T a {H{C F )) ^ 0. □ 



2 Fractal-Line Intersection 
2.1 Exact Intersection 

Theorem 2.1 A hyperplane intersects a hyperdense fractal if and only if it intersects its 
convex hull. This equivalence holds only if the fractal is hyperdense. 

Proof The proof is based on Cantor's Intersection Theorem and the address generation of 
F = (Ti, . . . ,T n ) C M. d in Theorem 1.2. We will show that the hyperdensity of F implies a 
decreasing sequence of compact sets, which tend to a point by Cantor's Intersection Theorem. 
The index sequence itself will correspond to an address, which in the limit locates a fractal 
point, since the fractal is the closure of all possible addresses. Let us now begin the proof. 
If a hyperplane L intersects F, it must clearly intersect Cp := Conv(F). On the other hand, 
if L intersects Cp, by the fractal's hyperdensity, it also intersects H(Cp). So L must intersect 
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T kl {Cp) for some ki G Af. Let this intersection be denoted as 1\ := Lfl T kl {Cp) C M. d . Then 
Ji is compact, because Cp is compact. 

Since T kl (H(Cp)) = lX=i °T k (Cp), according to Lemma 1.1 with a = the fact that 
L intersects T kl {Cp) implies that it also intersects T kl (H{Cp)) and thus T kl o T k2 (Cp) for 
some k 2 G Af. Once again denoting I 2 := L HT kl oT k2 (Cp) we have that this set is compact, 
and I 2 C LC]T kl (C F ) = h- 

Continuing to apply Lemma 1.1 in the above recursive procedure, by induction we get a 
strictly monotonically decreasing sequence (since T k , k G Af are contractive) of compact 
sets in K^: Ii D I2 D . . . D Ij D . . . each with a corresponding address composition: 
T kl o . . . o T k . . According to Cantor's Intersection Theorem D^Li h 7^ ^ an< ^ ^ contains a 
single point / G M d , since the address composition contracts to a point in the limit. 
Starting with any p G V HCf we have that / = lim^oo T kl o . . .0 T k , (p) G F by the address 
generation of F, so L intersects F in /. Note that such an intersection may not be unique, 
since in our recursive proof, we only chose one indexed set in the Hutchinson union at each 
step, though L may intersect multiple. 

The above is under the condition that F is hyperdense. If it is not, then by definition there 
is a hyperplane which intersects Cp but not H(Cp)- Since F C H(Cp) this hyperplane will 
not intersect F, countering the equivalence. □ 

Corollary 2.1 A line intersects a chain fractal in R 2 if and only if it intersects its convex 
hull. 

Algorithm 2.1 (Fractal-Line Intersection) Let F = (Ti, . . . ,T n ) be a hyperdense fractal 
(possibly a chain fractal), and assume that its convex hull Cp is known explicitly. Further- 
more, let L be our line, and e > an arbitrary stopping parameter. 
Step 0: a := (0), I := 
Step l:IfLn T a {C F ) ^ then 
For k = 1 to n 

If Lf]T a oT k (C F ) ^ then 
a := (a, k) 

7/diam(L n T a (C F )) > e then 

Go To Step 1 
Else 

Go To Step 2 
End 
End 
End 
End 

Step 2: I := IU(LnT a (C F )) 

Proof It follows directly from the proof of Theorem 2.1. The algorithm finds a union of 
intervals - each with length below e - that covers the points of intersection of the line and 
the fractal. □ 
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2.2 Approximate Intersection 



We now turn to examining the second question posed earlier, namely the number of inter- 
sections within some e > accuracy. These results are also relevant for applied reasons, 
such as the signal-reception efficiency of fractal antennas. If an antenna is designed so that 
it is hyperdense - possibly a chain fractal - then not only will it intercept all signal planes 
crossing its convex hull - making it space-efficient - but it will do so an infinite number of 
times, in an approximate sense. 

Theorem 2.2 If an open set has a common point with an IFS fractal, then it has an un- 
countably infinite number. 

Proof Let S C C the open set and / G FnS. We may suppose that / has a finite address 
a G A, \a\ < oo, since in the address generation of F, the fractal points with a finite address 
are dense in F, which is their closure according to Theorem 1.2. So we may replace / with 
another fractal point /' that has a finite address, and which is close enough to / to be still 
be an element of 5". Thus we may suppose that / has a finite address, with a corresponding 
map T a and primary fixed point p, meaning / = T a (p), p G V. Denote the contraction 
belonging to p as T G T. 

Let e' > be the radius of some ball centered at / - denoted by B' := B(f,e') - that is still 
contained in S. If we show that there are an uncountably infinite number of fractal points 
in B', then that implies the theorem. 

Let us map back / = T a (p) and B' by T~ x to p = T~ 1 (f) and B" := T~ 1 (B') = B(p,e") 
respectively. Here e" = where (p a is the product of the factors of the contractions in T 
making up T a . Then since F is compact, we may map it iteratively by T until it is contained 
in B", that is 3k G N : T k (F) C B", and clearly T k (F) has an uncountably infinite number 
of points. 

Mapping it all back by T a , we have that T a o T k (F) (ZFnB'cFnS. Thus we have shown 
an uncountably infinite number of common points of F and S. □ 

Corollary 2.2 Suppose that some open e > translational neighborhood L\ of a line L 
contains a point in an IFS fractal F, meaning 

Be > 0, / G F : / G L\ := {z G C : d(L, z) < e} 

Then there are an uncountably infinite number of fractal points in L\. 

Corollary 2.3 Suppose that some open e > angular neighborhood L" of a line L around 
q G L contains a point in an IFS fractal F , meaning 

Be > 0, / G F : / G L a £ := {z G C : Z(L, z - q) < e} 

Then there are an uncountably infinite number of fractal points in L a z . 

Both of the above corollaries follow directly from the theorem. The next theorem is of a 
slightly different nature, and is not implied by it. We restrict ourselves to IFS fractals in the 
complex plane, generated by an IFS of similitudes. 



7 



Theorem 2.3 Let T be an IFS with at least one map T G T having a rotation angle ft for 
which y~ is irrational. Then any line L that intersects F = (T) in some point f with a finite 
address, intersects it in at least a countably infinite number of points, with any e > angular 
accuracy around f , meaning 

^ > 3(/ fe )- =1 C F : Z(LJ k - f) < e, k G N 

Proof The address of / is finite, meaning that / = T a (p), p G V ', a G A, \a\ < oo. 
Let us transform back / and L to p G V by T^ 1 . Then we have that p = T~ 1 (f), and 
denote V := T~ l (L) which is also a line. The e angular neighborhood of L around / is 
transformed by T" 1 to an e angular neighborhood of V around p. Clearly p G V since 
/ = T a (p) G L = T a (L'). For the theorem to hold, it is sufficient to find an infinite number 
of fractal points within the angular neighborhood of L', since we can map these points with 
T a to the angular neighborhood of L. 

Let us choose any q G F, r G L' , with a := arg(g— p), (3 := arg(r— p). Then arg(T fc (g)— p) = 
(kft + a) mod 2ti and the iterations of q by T will be along a logarithmic spiral around p. 
We will show that the iterations visit the e angular neighborhood of V infinitely often. 




Figure 1: Illustration of the angular neighborhood of L' around p. 

We have supposed that ^ is irrational, so by the Equidistribution Theorem, the sequence 
{kft mod 27r)fc £ N is uniformly distributed in [0, 2tt), and thus it is also dense in this interval. 
So approximating the angle between q and r with respect to p, ie. (3 — a, we have the 
following: 

Vj G N 3kj G N : {{kjft + a) - p) mod 2vr < - < e 

Therefore the sequence (T kj (q))jen will be within the e angular neighborhood of V with 
respect to p. Thus mapping the sequence back by T a , it will be in the required e angular 
neighborhood of L with respect to /, and clearly (T a o T kj (q))j € fq C F. So with fj : = 
T a o T fcj (g), j G N we have shown a countably infinite number of fractal points in the e 
angular neighborhood of the intersecting line. □ 
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Corollary 2.4 Suppose a line intersects a planar IF S fractal in a point with a finite address, 
and has some e > angular neighborhood around the point that contains only a finite number 
of fractal points. Then all contractions in the IFS have a rotation angle •§ for which ^ G Q. 

The above corollary hints at the relevance of a certain class of planar IFS fractals, having 
roots of unity as rotation factors. They may hold a special place in a theory of the distribution 
/ connectivity of fractal sets. Examining such "rational fractals" further may prove to be 
a fruitful venture, and the case of primitive roots of unity may be even more worthwhile. 
Indeed these theorems seem to call for an investigation into the translational and angular 
distribution of IFS fractals. 

2.3 Quantifying Intersection 

In the previous section we have shown an infinite number of approximate intersections be- 
tween a fractal and a line under certain conditions. This implies the problem of determining 
how much of the fractal falls within a translational or angular neighborhood of the line. 
Clearly what we need is a probability measure which can be easily evaluated in practice. 
We will discuss the well-known IFS invariant measure, and provide a new formula for its 
computation. 

Definition 2.1 We say that a measure fi : M. d — > [0, 1] is invariant with respect to the IFS 
T = {Xi, . . . , T n } if it satisfies the following equation 

H(S) = w lf i(Tr 1 (S)) + ...+ w nf i(T-\S)) 

for any S C M. d in its domain, with some fixed weights Wk G [0, 1], Y2k w k = 1- 

When the IFS maps are similitudes of the form Tk(z) := pk + XkRk(z — pk) with G (0, 1) 
and unitary Rk G M dxd , then there exists a unique s > 0, called the similarity dimension of 
the fractal, for which Y2k ^1 = 1- ^ n suc h a case, the weights Wk = Aj? seem natural. 

Theorem 2.4 (Hutchinson [13]) With respect to any IFS and weights, there exists a unique 
invariant probability measure with bounded support. We call this the invariant measure with 
respect to the IFS. For any initial probability measure with bounded support fio, the recursion 
f^L = uii/jl-1 ° T{~ + . . . + WnfiL-i ° T~ x , L G N tends to the invariant measure as L — >■ oo. 

Theorem 2.5 For any IFS T = {Ti, . . . ,T n }, with weights w±, . . . ,w n , and primary fixed 
point p G V, the invariant measure v has the form 

v{S) = lim V(> a : a G A, \a\ = L,T a (p) G S) (S C R d ) where w a : = 

k=l 



W, 



a(k) 
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Proof The functions u := 5 P and Vl{S) '■= J2( w a '■ a E A,\a\ = L, T a (p) E S), L E N, S C 
M. d clearly satisfy the required properties of a measure, considering that 

1 = (wx + . . . + w n ) L = ^2(w a : a E A, \a\ = L) (L E N) 

The Dirac measure v§ has the support {p} which is bounded. We now show that the weighted 
recursion in Theorem 2.4 holds for ul, L E N. The key to its derivation, is breaking 
up the sum according to the first coordinate of each address, which is taken from M = 
{l,...,n}. 

Vl{S) = ^(wa : a E A, \a\ = L,T a (p) E S) = 

n 

= Yl Z)K a — (k,b) E J\f x A,\b\ — L — l,T k o T b (p) E S) = 

k=l 

n n 

= Y, Wk J2( w » --bEA,\b\=L- l,T b (p) E T k \S)) = J2 w ^l-i(T^(S)) 

k=l k=l 

Note that the derivation also holds for L = 1. By the above recursion and Theorem 2.4 we 
have that 31im(z/i) =: v and it is the invariant measure with respect to the IFS T . □ 

The above explicit formula for the invariant measure has clear theoretical and practical 
advantages over existing methods for its computation. One method well-known in image 
processing, is to discretize the plane over a rectangle, and approximate the invariant measure 
via a matrix recursion. Another method is carried out stochastically using Elton's Ergodic 
Theorem [5] . Our exact formula for the invariant measure v assigns weights w a to each fractal 
point T a (p) with an address of length L, and checks which points fall into the set S, then sums 
the weight of those points. Clearly v can also be approximated in practice via its recursion 
up to Ul, with a large enough level L E N. We may thus find the density of intersections 
with the fractal in between parallel rays of light, reasonably spaced at e := A^diau^'P), to 
see the intensity of the shadow in a particular direction, as illustrated on the figure below. 
This kind of analysis has clear applications to the design of fractal antennas for instance, 
and the tomography of fractalline structures. 




Figure 2: Ray absorption density plots under rotational perturbation. 
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3 Conclusive Remarks 



In our paper, we have highlighted the connection between the projection of IFS fractals and 
the Fractal-Line Intersection Problem, relevant to applications in computer graphics and 
antenna design. We have introduced classes of fractals for which the projection from any 
direction is a segment in the plane. We have done this in a general setting, examining the 
intersection of hyperplanes and IFS fractals in R d , del Our method implies an algorithm 
for finding the points of intersection. We have discussed the cardinality of the intersection, 
and gave an explicit formula for the well-known invariant measure of IFS, in order to quantify 
intersection. 

We hereby express our gratitude for the support of the Faculty of Mathematics at the Uni- 
versity of Waterloo, and primarily that of E.R. Vrscay who is supported by NSERC. 
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